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ABSTRACT 

We apply the optimal filter technique to Sloan Digital Sky Survey photometry around 
Segue 1 and find that the outer parts of the cluster are distorted. There is strong ev- 
idence for ~ 1° elongations of extra-tidal stars, extending both eastwards and south- 
westwards of the cluster. The extensions have similar differential Hess diagrams to 
Segue 1 and a Kolmogorov-Smirnov test suggests a high probability that both come 
from the same parent distribution. The location of Segue 1 is close to crossings of 
the tidal wraps of the Sagittarius stream. By extracting blue horizontal branch stars 
from Sloan's spectral database, two kinematic features are isolated and identified with 
different wraps of the Sagittarius stream. We show that Segue 1 is moving with a 
velocity that is close to one of the wraps. At this location, we estimate that there are 
enough Sagittarius stars, indistinguishable from Segue 1 stars, to inflate the velocity 
dispersion and hence the mass-to-light ratio. All the available evidence is consistent 
with the interpretation that Segue 1 is a star cluster, originally from the Sagittarius 
galaxy, and now dissolving in the Milky Way. 
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1 INTRODUCTION 



Segue 1 was discovered by IBelokurov et al.l l|2007f ) as an 
overdensity of resolved stars in imaging data from the Sloan 
Digital Sky Survey (SDSS). It is located at equatorial co- 
ordinates 02000 ~ 152°, ^2000 ~ 16° and has a heliocentric 
distance of 23 ± 2 kpc. This corresponds to Galactic coor- 
dinates i — 220.5°, b = 50.4° and a Galactocentric distance 
of ~ 28 kpc. IBelokurov et al.l (|2007h suggested that Segue 1 
was an extended globular cluster, possibly associated with 
the Sagittarius stream. With its initially determined half- 
light radius of approximately 30 pc, it would be amongst 
the largest Milky Way globular clusters such as Palomar 5. 
The object is also unusually faint {My ~ —3) for its size. 

Segue 1 has some points in common with two 
ot her discoveries mad e with SDSS - namely Willma n 
1 (| Willman et al.ll2005D and Bootes II l| Walsh et ahlbOOTh . 
All three objects have similar absolute magnitudes and half- 
light radii, intermediate between the dwarf spheroidal galax- 
ies, which are dark matter dominated, and the globular clus- 
ters, which show no evidence for dark matter (see Figure 8 in 
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l| Walsh et al.ll2008h ). Determining the true nature of these 
three objects is important as it may shed light on the impor- 
tant question of the size of the smallest dark matter haloes 
in which baryons collapsed to form galaxies. 

The interpretation of Segue 1 as a glo b ular c luster 
has recently been contested by iGeha et al.l (|2009| ). Us- 
ing Keck/DEIMOS spectroscopy, they measured the ra- 
dial velocities of 24 stars in Segue 1 with a mean helio- 
centric velocity of ~ 206 kms~^ and a velocity dispersion 
of 4.2 ± 1.2 kms"'^, leading to claims that Segue 1 is a 
dwarf galaxy rather than a globular cluster. Assuming these 
stars are gravitationally bound to Segue 1 and are in dy- 
namical eq uilib rium, then the implied mass-to-light ratio 
is ~ 1200 (|Geha et all 120091 ') . which would make Segue 1 
the most dark matter dominated galaxy detected to date. 
However, both assumptions are questionable. If Segue 1 is 
a globular cluster that is undergoing tidal disruption, then 
extra-tidal stars may not be so easy to distinguish from grav- 
itationally bound members and the hypothesis of dynamical 
equilibrium may be a poor one. More seriously, if Segue 1 is 
immersed in the Sagittarius stream, then contamination of 
any sample of Segue 1 stars by stream stars may be hard to 
avoid. 
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In this paper, we analyze SDSS and auxiliary Canada- 
France-Hawaii Telescope (CFHT) photometry of Segue 
1. We use the optima l filter techniques pioneered by 
lOdenkirchen et al.l jioosh to identify tidal features in §2 and 
§3. A stellar population embedded deep in a massive halo 
would not be expected to show visible signs of battering by 
the Galactic tides. We compute Segue I's structural param- 
eters in §4. Finally, in §5, we extract from the SDSS spectral 
database blue horizontal branch stars and used them to iden- 
tify the kinematical signal from the Sagittarius stream. We 
show that contamination of kinematically selected Segue 1 
stars by Sagittarius stream stars is a serious problem, and 
can lead to artificially inflated velocity dispersions. 



2 DATA 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) is 
an imaging and spectroscopic survey that covers one-quarter 
of the celestial sphere. The SDSS data are described in the 
data release papers (Adelmann- McCarthy 2008 for the sixth 



release, DR 6) and documented at ' [http:/ /www.sdss.org ' 

We select stars in a 10° x 10° box centered on Segue 
1 from SDSS DR 6 data, using the SDSS clean photom- 
etry flags and a magnitude cut (r > 22) to remove data 
reduction artifacts from the fi eld. In addition, we correct for 
extinc tion using the maps of lSchlegel. Finkbeiner. fc Davis! 
dlOQSt) and apply the Ube rCal correction as described in 
[Padm anabh an et al.l |20o3). Figure [1^ shows the density 
of all stars (fainter than r — 14) selected by dividing the 
area in 75 x 75 pixels and smoothing the field with a Gaus- 
sian kernel FWHM of 2 pixels (this is applied to all den- 
sity plots in this work). Segue 1 is too faint to stand out 
in this figure. Figure [TJd shows a density plot of stars that 
we have removed from the sample using our magnitude cut. 
Structural artifacts due to problems in data reduction are 
clearly visible and correspond to the SDSS scan patterns 
on the sky. Figure [ij; shows our estimate of the field star 
density with extinction contours overplotted. The estimate 
is determined by first replacing Segue 1 within 1° with a 
representative patch of the field star distribution taken at 
Q = 149° , S = 19° . By a similar cloning method, we remove 
the other obvious overdensity in the field, the Leo I dwarf 
galaxy {a ~ 152.2°, (5 ~ 12.5°, at a Galactocentric distance 
of 250 kpc). We then compute the density and smooth the 
resulting distribution with a Gaussian kernel with a FWHM 
of 5 pixels and a box-car smoothing over 25 pixels. Finally, 
since classification of galaxies and stars especially at faint 
magnitudes can be uncertain in the SDSS, it is important 
to verify that possible structures are not influenced by such 
misclassifications. Figure [TJi shows that there are no obvi- 
ous galaxy overdensities in the immediate vicinity of Segue 
1. There does seem to be an underdensity, however, this 
would not generate spurious tidal structures near Segue 1 
except in the highly unlikely event that almost all galaxies 
were misclassified as stars. 

The main dataset used in this paper is from the SDSS. 
However, we also have access to a set of MegaCam pointings 
taken on the Canada-France-Hawaii Telescope (CFHT) on 
17-25th January 2007. Data were taken in g and i filters, 
typically three 200 s exposures using the measured zero- 
points for MegaCam. The MegaCam g and i magnitudes 



Table 1. K-S tests comparing the Hess diagrams of different 
above average signal regions and Segue 1. A high K-S probability 
indicates a high confidence that the two samples are drawn from 
the same parent population. The number of stars is the count of 
stars inside the CMD mask in a given box. The count of Sagit- 
tarius stars (i.e. stars falling inside the CMD mask but belonging 
to the Sagittarius stream rather than Segue 1) is estimated using 
the methods described in 55.1. 



Debris 


Area 


Number 


Number 


K-S 


Region 


[deg2] 


of Stars 


of Sag Stars 


Probability 


Segue 1 


0.045 


59 


11 - 25 




Box 1 


1.140 


1175 


285 - 627 


0.93 


Box 2 


1.785 


1908 


446 - 982 


0.99 


Box 3 


0.515 


607 


128 - 283 


0.99 


Box 4 


0.825 


833 


206 - 453 


0.06 


Box 5 


0.350 


370 


87 - 193 


0.55 


Box 6 


0.600 


592 


150 - 330 


0.15 


Box 7 


0.440 


391 


110 - 242 


0.05 



are converted to SDSS magnitudes by matching the stars in 
the CFHT data to the SDSS observations. We find 21952 
matches between the two data sets. We then plot the differ- 
ence in magnitudes in the two systems versus CFHT color 
to find the relations 



5SDSS — ffCFHT 



*SDSS — iCFHT 



-0.22 + 0.21(gcFHT - icFHT) 

— 0.05((;cFHT — jcfht)^ 
-0.09 -I- 0.01(gcFHT - icfht) 

+0.01((/CFHT — JCFHt)^ 



(1) 



There are 10 MegaCam fields each covering 1 deg , strad- 
dling Segue 1 at fixed declination. A rectangle bounding the 
10 fields is shown in Figure [T^. Although deeper than the 
SDSS data by about a magnitude, the CFHT data cover 
considerably less area and are not so useful for diagnosing 
tidal features. Nonetheless, it is the best dataset for estimat- 
ing the surface brightness profile and hence the structural 
parameters of Segue 1. 



3 THE OPTIMAL FILTER TECHNIQUE 
3.1 Introduction 

The optimal filter technique works by calculating condi- 
tional probabilities of satellite and foreground membership 
from densities in colour-magnitude space, known as Hess di- 
agrams. For each star, the likelihood (or weight) of Segue 1 
membership is simply proportional to the ratio of the satel- 
lite and field Hess diagrams in the relevant pixel of colour- 
magnitude space. As the Galactic foreground contains stars 
from all possible populations in a large range of distances, 
the membership likelihood never reaches certainty. There- 
fore, instead of looking at individual stars, we study the den- 
sity distribution of possible members represented by weights 
summed in pixels o n the celes tial sphere. The method is de- 
scribed in Odenkirchen et al.l (j200^) and summarized by the 
formula 



nc(fc) 



S, [n(fc, j)/c(j)//F(i) - np{k,j)fcU)/fF{j)] 



(2) 



S./S(j)//f(j) 
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(a) (b) (c) (d) 

Figure 1. Panel (a) shows the distribution of all stars with clean SDSS photometry with 14 < r < 22 in a 10° X 10° square centered 
on Segue 1 (dark areas are high density). The overdensity at a si 152.2°, (5 12.5° is the Leo I dwarf galaxy. Also shown as a black 
rectangle is the CFHT footprint, while the white circle marks the location of Segue 1. Panel (b) shows the density of stars which have 
been removed from the sample using our magnitude cut. The visible bands correspond to the SDSS scan directions. Panel (c) shows our 
estimate of the field stars (background and foreground) with extinction contours overplotted. The sidebar gives the number of stars per 
square degree. As will become clear later, there does not seem to be any correlation between high extinction and any extra-tidal features. 
Panel (d) shows the distribution of galaxies from the SDSS data in our field of view. The sidebar shows number of galaxies per square 
degree. 




g-r g-r 



(a) (b) (c) 

Figure 2. Panel (a) shows a scatter plot of stars within 0.12° of Segue I's center, as seen in the CFHT data. The M92 ridgeline is 
plotted and it is in very good agreement with the data, which suggests that we use the more complete data on M92 to approximate 
the CMD of Segue 1. The mask overplotted is used to select possible members. Stars outside the mask are highly unlikely to be cluster 
members and their weight is zero. Panel (b) shows the Hess diagram of M92. A spread has been added to the data in order to more 
realistically reproduce the SDSS photometry. The spread corresponds to the difference between the uncertainty inherent in the M92 data 
and the uncertainty found in SDSS at the corresponding magnitudes. Panel (c) shows the ratio of the M92 Hess diagram and the field 
Hess diagram. This ratio is used as the weights distribution in the optimal filter technique. 



where nc{k) is the background and foreground corrected, 
weighted Segue 1 density in position space; n{k,j) labels 
the stars that are in the j*'' bin in color-magnitude space 
and in the fe*'' bin in position space; fc and fp represent 
the color-magnitude Hess diagrams of Segue 1 and the field 
respectively (their ratio being the weight we employ). Fi- 
nally, np{k,j) is the field star density, which we show in 
Figure [!};, split according to the bins in colour-magnitude 
and position space. 

We analyzed the color-magnitude distribution of Segue 
1 (fc) by considering all stars that lie within a 0.12 deg aper- 
ture around the cluster center. Th e Hess diagram o f Segu e 
1 is shown in Figures 3 and 4 of Belokurov et al.l (|2007l ). 
with the ridgeline of M92 from lClemI (j2005l ) at the distance 
modulus 16.8 mag overplotted. We determine the field star 
color-magnitude distribution (fp) by considering stars out- 



side of an aperture of 0.4°. We find that optimal results are 
achieved by subdividing color-magnitude space with 50 by 
50 pixels and smoothing the resulting distribution with a 
Gaussian kernel with a FWHM of 3 pixels in both cases. In 
the optimal filter technique, the ratio of the two Hess dia- 
grams is used as the weight to determine whether a star is a 
member of Segue 1. Guided by the ridgeline, a mask about 
the relevant region of the Hess diagram is tightly drawn. 
The weight distribution which we derive using this method 
is rather choppy, due to the small number of stars found close 
to the clusters center (< 100 stars) and such a distribution 
seems unphysical. 

Using the deeper CFHT data, we find that the agree- 
ment between the Segue 1 CMD and the M92 ridgeline is 
sufficiently close (Figure [2] a) to approximate the CMD of 
Segue 1 with that of M92. Figure [2d shows the M92 Hess 
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Figure 3. Panel (a) shows the weighted and background corrected Segue 1 star density with a structure extending outwards from Segue 
1. The density is normahzed to the average pixel value within 0.1° of the center of Segue 1. The contours show 1.5,2,3,4 and 5 a levels 
above average density, whilst the colour within the contours can be converted into normalized pixel density by the key. Panel (b) shows 
the distribution of pixel values for panel (a) with a fitted Gaussian. The vertical colour bars indicate the 1.5,2, 3, 4, and 5fT levels as 
determined from the Gaussian distribution. The Gaussian fits the distribution well, but there is an excess at the positive end due to the 
presence of Segue 1. 



diagram, with a spread added in order to more accurately 
reflect the SDSS photometry, as well as the ratio of M92 
and field Hess diagrams which we ultimately use as weights 
(Figure [2};). It is noticeably smoother than the weights dis- 
tribution generated by using SDSS data alone. 

The foreground is composed of the smooth halo popu- 
lation and some Sagittarius material. While the mean chem- 
ical p roperties of the oute r halo are not very different from 
M92 (|Carollo et al.ll2007D . halo stars are of course much 
more spread out both spatially and in metallicity. Sagit- 
tarius stars at this location occupy a slightly brighter and 
redder part of the CMD, because Sagittarius is more metal- 
rich. Nonetheless, there is an overlap between Sagittarius 
stars and Segue 1 stars on the CMD. In other words, an 
M92 mask centered on Segue I's sequence will pick up some 
Sagittarius stars. 

3.2 Segue 1 as seen by the Optimal Filter 

Using the optimal filter technique with the stellar weights 
described above, we find an extended structure which ap- 
pears to be connected to Segue 1. In Figure [3^, we show 
the stellar density nc- The filled coloured contours indicate 
the 1.5,2,3,4, and 5 a levels above average density. There is 
a clear structure extending in a horizontal band from the 
center of Segue 1 out to approximately a = 154°, 5 = 15.8°. 
It is approximately 0.5° wide, and lies along the direction of 
the Sagittarius tidal tail at this location. In addition, there 
appears to be a structure extending downward from the clus- 
ter. Figure|3)D shows the distribution of stellar weight counts 
in the pixels. We associate the excess of the distribution 
above the fitted Gaussian at the positive end with Segue 1 
and its tidal debris. 

We repeat the optimal filter technique defining the stel- 



lar weights using the g — i versus i Hess diagram as well 
as the ci vers us i Hess diagram. Here, ci is a color index 
introduced by lOdenkirchen et al.l l|2002h and defined by us 
as: 

ci = 0.918(3 -r)+0.397(r- i) (3) 

This color index is chosen to lie along the one-dimensional 
distribution of M92 stars in the g — r versus r — i space. We 
note that the stars which we select from the SDSS as being 
in the center of Segue 1 do not lie along a one-dimensional 
locus. We attribute this to contamination from foreground 
stars. By creating this color index, the greatest amount of 
available data is used and hence the results of the optimal 
filter technique based on it should be most robust. In Figure 
3] we show 4 different views of the optimal filter technique 
with weights derived from g ~ r versus r, g — i versus i and 
ci versus i Hess diagrams. In addition, we show the average 
of the g — r versus r and g — i versus i based optimal filter 
results. 

3.3 Extra-Tidal Features 

It is essential to verify that the visually identified structures 
are actually related to the cluster and not merely the result 
of chance alignments with noise within the field. To this, end 
we analyze the Hess diagrams of a number of overdensities 
in the field. 

We investigate 7 areas of overdensity near Segue 1 in 
more detail. For each of the boxes shown in Figure |31 we 
generate differential Hess diagrams in order to determine if 
the stellar populations are similar to that of the core of Segue 
1. The lower panels of Figure 3] show some sample results. 
We find that boxes 2 and 3 (which represent tidal debris 
directly connected to the cluster) have differential Hess di- 



© 2009 RAS, MNRAS 000.[HfT2l 



The Origin of Segue 1 5 




J I I I I I I I I I I I I I J I I I I L 

156 154 152 150 148 

RA [Deg.] 




-0.5 0.0 0.5 1.0 1.5 -0.5 0.0 0.5 1.0 1.5 -0.5 0.0 0.5 1.0 1.5 -0.5 0.0 0.5 1.0 1.5 -0.5 0.0 0.5 1.0 1.5 
g— r g-r g-r g— r g-r 



(e) 

Figure 4. Panel (a) shows the result of the optimal filter technique with weights defined from the g—r versus r Hess diagram. Panel (b) 
shows the results using weights from g—i versus i. Panel (c) is the average of the figures shown in panels (a) and (b). Panel (d) shows the 
results of the optimal filter technique using weights from the ci versus i Hess diagram. Although there are slight differences between our 
different implementations of the optimal filter technique, all four plots support that there are extra-tidal features around Segue 1. Panel 
(e) shows differential Hess diagrams of Segue 1 and boxes 2-5 respectively. Boxes 2 and 3 are very similar to Segue 1 whereas boxes 4 
and 5 are not, even though they have reasonably high significance in the optimal filter analysis. 
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agrams that follow the theoretical ridgeline of M92 closely. 
The slight deviation from the ridgeline at faint magnitudes 
is attributed to deficiencies in the SDSS photometry (also 
visible for example in the Hess diagram of Belokurov et al. 
2007). Boxes 4 and 5 lie close to the cluster and have a high 
significance. However, we find that their differential Hess di- 
agrams do not agree with that of Segue 1. Box 1 is a high 
signal area that lies on what seems to be an extended row 
of overdensities near Segue 1, but again seems to have a dis- 
similar differential Hess diagram. Boxes 6 and 7 represent 
two overdensity patches clearly below the cluster (with box 
7 only appearing as a significant overdensity in the g — r and 
g — i implementations of the optimal filter technique) . Their 
differential Hess diagrams do not look like those of Segue 1. 

We assess the similarity between the differential Hess 
diagrams of the boxes with that of Segue 1 by using the 
Kolmogorov-Smirnov (K-S) test. In order to apply the K-S 
test, we in effect determine a luminosity function from the 
g — r versus r Hess diagrams by summing pixel rows and 
combining the counts in adjacent rows (i.e. doubling the 
size of the pixels in the magnitude direction). To reduce the 
influence of foreground stars on our test, we exclude noise 
seen in the Hess diagrams red ward of g — r = 0.8. Using 
only stars in a narrowly defined mask about the ridgeline 
proved impractical because the noise in the Hess diagrams 
is too great. We calculate cumulative luminosity functions 
and determine the maximum difference between the cumu- 
lative distribution found in the box to that of Segue 1. The 
results of the K-S tests are summarized in Table [T] with a 
high probability suggesting that the two distributions come 
from the same parent distribution. The table also gives the 
number of stars within each box, together with an estimate 
of the number of Sagittarius stream stars (see Section 5.1). 
We caution that smoothing the data may have introduced 
correlations and so the K-S probabilities may be slightly too 
large. We conclude that the overdensities in boxes 2 and 3 
are consistent with being tidal material stripped from Segue 
1, at least as far as the Hess diagrams are concerned. 

Finally, we also investigate whether or not the extra- 
tidal features visible in Figure |4] are the result of chance 
alignments with the noise in our field of view. To this end, 
we place 4424 simulated Segue 1 like objects in the field. 
These are created by randomly sampling the Plummer pro- 
file fit to the CFHT data and the luminosity function gener- 
ated from M92. We perform the optimal filter analysis and 
measure the stellar excess around Segue 1 and the simulated 
object by counting the stellar weights in a circular annulus 
between 0.3 and 1 degree from the cluster center, as indi- 
cated in Figure [S] We count only those pixels that are at 
least 1.5(T above the average density. The excess measured 
around Segue 1 is a significant outlier (at least 3(t) com- 
pared to the excess measured about the simulated objects, 
as shown in Figure |6l 

Note that the noisy background in Figures |3] and 2] 
seems to be a consequence of the location of Segue 1 in the 
Sagittarius stream. Suppose we place a simulated Segue 1 in 
a patch of the sky at equatorial coordinates a = 152° , 5 = 
46° , which is well away from the Stream. Applying the opti- 
mal filter technique yields the results shown in Figure [5] The 
object is cleanly and easily recovered and the background 
has very little prominent substructure. This suggests that 
the "noise" in our earlier Figures is probably real substruc- 




156 154 152 150 148 
RA[Deg.] 

Figure 5. The optimal filter technique is applied to input data in 
which a simulated Segue 1 lies in a part of the sky uncontaminated 
by the Sagittarius stream. The field star distribution is much 
smoother and the noise is significantly reduced as compared to 
Figures [3] and m 



ture in the Sagittarius stream, although further investiga- 
tions are need to confirm this. 

To conclude, we have demonstrated that there are struc- 
tures (e.g., boxes 2 and 3) seemingly connected to Segue 1 
and which have very similar differential Hess diagrams to 
Segue 1. The structures are unlikely to be caused by chance 
alignment with noise generated by the Sagittarius stream. 
It is very natural to ascribe them to tidal features. Note 
that Martin et al (2008) have argued that, even with the 
smoothing that tends to obliterate small-scale structure, re- 
alizations of purely spheroidal models with small numbers 
of stars can appear just as distorted as the observed pho- 
tometric data on some of the ultra-faint galaxies. However, 
their calculations do not directly address the existence of 
tidal features in the outer parts of these objects. 



4 STRUCTURAL PROPERTIES OF SEGUE 1 

Using available CFHT data, we are able to get a view of 
a smaller, but deeper, field containing Segue 1. Performing 
the optimal filter analysis on the data proved difficult since 
the field is too small to define a meaningful background. We 
are only able to identify the cluster but not any prominent 
extra-tidal features found using the SDSS data. We are able 
to confirm that this is not due to differences or anomalies 
in the SDSS and CFHT data by applying the optimal filter 
technique to a field the size of the CFHT field of view cut out 
of our SDSS data. Doing this, we find a similar result that 
the cluster can be isolated, but not any obvious extra-tidal 
features. 
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Figure 6. The distribution of excess measured around the sim- 
ulated Segue 1 objects placed randomly in the field of view. The 
horizontal axis show the excess as a normalized count i.e. the ra- 
tio of the count around the simulated object to the count around 
Segue 1. The vertical axis indicates how many simulations had 
this excess. The dashed line indicates the location of the excess 
around Segue 1. Depending on the Hess diagram used to define 
the weights in the optimal filter technique, the excess around 
Segue 1 is at least a 3(7 outlier. Here we show the results of the 
simulation with weights defined by the ci versus i Hess diagram. 



With the CFHT data, we analyze the cluster density 
profile and luminosity in greater detail. Figure [7] shows 
Plummer, King, and exponential profiles fit to the surface 
density of Segue 1. The parameters determined from the fits 
are summarized in Table [2] and compared to the values de- 
termined earlier by Belokurov et al. (2006) and Martin et 
al. (2008). In the fitting, only those stars that fall inside 
the CMD mask of Segue 1 shown in Figure [2] are consid- 
ered. There is evidence for extra-tidal populations in the 
succession of datapoints that lie above the assumed fit in 
the outer parts. Even though the deviations from the fit are 
always within 1 a, there are ~ 10 such datapoints. Within 
3 half-light radii, the fit to the models (King, Plummer or 
Exponential) is excellent with a of 1-1, 0.9 and 0.9 respec- 
tively. However, if the datapoints between 3 and 7 half-light 
radii alone are used, the values of the fit are 4.7, 3.6 and 
3.7 respectively, emphasising that the outer parts deviate 
from the smooth model. 

As noted by Martin et al. (2008), determining the lu- 
minosity for faint satellites with few observed member stars 
is challenging as it is strongly dependent on the inclusion 
or exclusion of a small number of stars that have evolved 
high up the red giant branch. We determine the luminos- 
ity of Segue 1 in three different ways, each using the CMD 
mask: (1) counting the flux within the Plummer half-light 
radius determined from our fits and doubling this to get the 
total luminosity, (2) determining the flux within the half- 



Radial star Count Distribution 




0.01 0.10 1.00 

Radius of Inner anulus in degrees 

Figure 7. The points show the number of stars in successive 
annuli about the cluster center (annuli width of 0.025°) with error 
bars assuming Poisson statistics. We only count stars that fall 
inside the CMD mask. The overplotted lines show King (solid 
line), Plummer (dashed line), and Exponential (dash-dotted line) 
fits to the distribution. Using the fits we determine a number of 
statistics for Segue 1 which are summarized in Table [2] Note the 
succession of points la away from the fits in the outer regions of 
Segue 1, which could confirm oxtra-tidal populations. 

light radius and scaling this up to 10", (3) counting the flux 
within 10". In each case, we subtract an estimate for the 
field star flux and add in a proxy for the missing flux that 
is lost due to our magnitude cut-off. The fleld star flux is 
determined by considering stars that fall within our CMD 
mask and that are in areas near the cluster but that our 
optimal filter technique shows to contain few cluster stars 
(two boxes defined from a — 152.2° to a = 156.3° and from 
S = 16.1° to S = 16.6° and a = 149.3° to a = 151.1° and 
from 5 = 16.1° to 5 = 16.6°). The missing flux is deter- 
mined from the amount of flux found below our cutoff in 
M92. However, as our magnitude cut-off in the CFHT data 
is i = 23, very little flux is missed (less than 1%). The de- 
termined values are summarized in Table [2] 



5 KINEMATICS OF SEGUE 1 
5.1 The Sagittarius Connection 

Segue 1 lies in a very busy area of the sky with multi- 
ple wraps of the Sagittarius stream as well as other debris 
streams, as shown in the "Field of Streams" (Belokurov et 
al. 2006). It has been known since the work of Ibata et al. 
(2001) that the Sagittarius tidal debris is wrapped around 
the Galaxy a rmmber of times. There are at least 4 wraps 
correspond ing to two streams, leading and trailing (see e.g., 
iFellhauer et al.il2006l ). A schematic plot is shown in Figure |8l 
in which the orbit of the Sagittarius dwarf is integrated for 
r^ 2 Gyr backwards and forwards in time. Segue 1 is located 
close to an intersection of the trailing and leading arms. 
The "Field of Streams" traces the young leading (A) and 
old trailing (B) wraps, which are closely matched in dis- 
tance around the North Galactic Cap. Simulations suggest 
the existence of further streams corresponding to old leading 
and young trailing wraps in the same field of view. So far, 
observational evidence for these wraps is sparse. 
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Table 2. Structural Parameters of Segue 1. Subscripts P, K and E 
refer to the Plummer, King and Exponential Fits, whilst all uncer- 
tainties are determined via bootstrapping. 




^ , , , , \ , , , , \ , , I 364 ± 147 

-50 50 

X [kpc] 



Figure 8. Schematic plot of possible orbits of the Sagittarius 
dwarf in the (x, z) plane integrated forwards (dark grey) and 
backwards (light gray) for 2 Gyr. The position of Segue 1 is 
close to crossings of the rosette orbit. The potential is defined 
by equations (l)-(3) in Fellhauer et al. (2006). The dashed line 
corresponds to an orbit starting at right ascension a = 283.7°, 
declination 5 = —30.5° and a heliocentric distance of 25 kpc (the 
location of the Sagittarius) with initial velocities Vr = 137 kms~^, 
/ic = —3.02 mas yr~^, /if = —1.49 mas yr~^. The solid line starts 
at the same location but has initial velocities Vr = 130 kms~^, 
/ia = —3.44 mas yr~^, /i^ = —1.32 mas yr~^. 
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Figure 9. Left: CMD showing the selection box used to pick BHB 
stars in the SDSS spectral database and the mask to pick turn-off 
stars for measuring the density of the Sagitta rius stream around 
Segue 1. Possible Segue 1 members selected bv lGeha etHI (|2009| ) 
are shown as black dots. Right: The upper solid line shows the 
density profile of the turn-off stars selected, with the dotted line 
showing a linear fit to the background and the residuals shown as 
lower dashed line. We also estimate the background by examining 
stars at the same Galactic latitude but with the sign of longitude 
reversed. This results in a lower background estimate and a higher 
density for the Sagittarius stream (upper dashed curve) . The filled 
circles mark the location of Segue 1. The vertical dashed lines 
mark the boundaries of the Sagittarius stream. 



" The definition of surface brightness in Martin et al. takes into con- 
sideration the ellipticity of the cluster, whereas we calculate an ef- 
fective surface brightness as the total flux divided by the area within 
the plummer half-light radius. Using our method and Ly and Ni, 
from Martin et al. one would find fiy = 28.6. 

^ Note that the three different values quoted for the absolute mag- 
nitude, central surface brightness and total luminosity correspond to 
the three different methods described in the main text. 

At the location of Segue 1, the wraps have similar dis- 
tances and are composed of similar stars. The only way to 
identify them is through their different kinematics. To this 
end, we interrogate the SDSS spectral database to pick out 
blue horizontal branch stars. Figure[9]shows the BHB branch 
of M92 with an enclosed box, which is used for selection in 
g — r versus r. To eliminate false positi ves, we simultane - 
ously impose the cut 0.9 < u — g < 1.5 (jSirko et al.ll2004l ). 
This box selects BHB stars at the approximate distance of 
Segue 1. 

Our rationale is to trace with BHBs a slice through the 
Sagittarius stream at constant right ascension (145° < a < 
155°). Any features that stand out in the same range of 
declination can be attributed to the Sagittarius stream. The 
right panel tells us what this range of declination is. The 
density of turn-off stars is shown as the uppermost solid 
line, with the turn-off stars selected via the mask shown in 
the left panel. The mask is wide enough to include all stars 
(22, excl uding 2 BHBs) in the spectroscopically confirmed 
sample of lGeha et al.l (|2009l ) . Although the mask is centered 
on the CMD of Segue 1, because of its width, it will also pick 
out stars in the Sagittarius stream. 

To make Sagittarius stand out more clearly, we subtract 
the Galactic foreground. The number density of Sagittarius 
stars depends on the model adopted for the foreground. In 
one method, we take advantage of the apparent linearity of 
the foreground (shown as a dotted line). Once subtracted, 
this gives the lower dashed profile, which peaks at ~ 300 
stars per deg'^. However, at low declinations, the fit could 
be in error due to a contrib u tion from the "Vir go Overden- 
sity" iBelokurov et al.l (|2006l ) ; lJuric et ah I (|2008l , see e.g.). As 
an alternative, we examine the density of stars at the same 
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Figure 10. Upper: The density of BHB stars in the plane of 
declination versus radial velocity (heliocentric for the left, Galac- 
tocentric for the right panel). The efficiency of the spectral 
follow-up is indicated by the horizontal bands with white rep- 
resenting 50% and black representing 10%. To guide the eye, 
the one-dimensional density profile of the Sagittarius stream 
(dashed curve from Fig (9)1 is shown connecting the two panels. 
The dashed lines mark the boundaries of the Sagittarius stream. 
Lower: Radial velocity distributions for stars in the declination 
range 10° 5 ^ 36° of the Sagittarius stream. The stream is 
most clearly identifiable in the right panel in which the velocities 
are Galactocentric. Note the two features at i'gsr = —105 (the 
A and B streams) and 130 kms~^ (the C and D streams). We 
also show as dotted lines the characteristic velocity distributions 
of the Galactic halo and thick disk to emphasise that the features 
cannot be ascribed to these components. Finally, the two black 
distributions are Gaussian fits to the stream velocities. 



Galactic latitude but with the sign of longitude reversed. 
This results in a lower foreground estimate and a higher 
peak density for the Sagittarius of ~ 600 stars per deg~^. 
The Sagittarius stream is limited to declinations satisfying 
10° < 5 < 36°. 

Figure [10] shows the distribution of BHB stars in the 
plane of radial velocity versus declination. Once the gra- 
dients due to the projection of the Local Standard of Rest 
have been removed, the two overdensities coincident with the 
Sagittarius stream in declination stand out at vqsr ~ —105 
kms~^ and 130 kms~^. At this location, it is unlikely that 
these features can be due to the Galactic thick disk or halo, 
as evidenced by the velocity histograms in the lower panel. 
In fact, they nicely match the velocities of the wraps of the 
Sagittarius stream at this location shown in Figure 3 of Fell- 
hauer et al. (2006). Segue 1 is in the Sagittarius stream 
and appears to be moving with the same velocity as one its 
wraps. In the lower panels in Figure 1101 we show Gaussian 
fits to each of the overdensities associated with Sagittarius. 
The dispersions are ~ 20 kms~^. Assuming typical radial 
velocity errors of ~ 15 kms~^ gives the intrinsic dispersion 
in each wrap to be ~ 13 kms~^. Of course, this is approxi- 
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Figure 11. The measured velocity dispersion of Segue 1, as a 
function of the number of contaminating stars from the Sagit- 
tarius stream. The lines indicate the average of 10000 simula- 
tions, whilst the error bars represent the spread. The black line 
represents a contaminating population with dispersion 10 kms~^ 
whilst the grey line represents a contaminating population with 
dispersion of 15 kms"'^. The value reported by Geha et al. (2008) 
is shown as a horizontal dotted line, with the dash-dotted lines 
giving the reported error. 



mate, as both measured dispersion and the errors could be 
somewhat larger. 

This has immediate consequences for estimates of the 
velocity dispersion of Segue 1, as the problem of contam- 
ination by Sagittarius stream stars is substantial. Let us 
es timate the numbe r of Sagittarius stars in the sample 
of lGeha et al] (|2009l ). As shown in Figure H our mask al- 
ready wraps around the candidate stars of Geha et al. Hence, 
the total number of Sagittarius stars at the location of 
Segue 1 can be read off the dashed curve in the right-hand 
panel. It is between ~ 250 and ~ 500 stars per deg^ depend- 
ing on the model adopted for foreground subtraction. To es- 
timate how many of these move with vgsr = 130 kms~^, 
we note that the two velocity peaks in the lower panel of 
Figure [10] integrate to roughly the same numbers of stars. 
So, allowing a factor of 0.5 to to account for the fact that 
both leading and trailing arm stars occur at this location 
and scaling to the 0.03 deg^ field of view used by Geha et 
al, we estimate that the final contamination in their sam- 
ple is at least ~ 250 stars deg~^ xO.5 x .03 deg^, which 
comes to 4 stars. This might be raised as high as 8 stars, if 
the lower value for the foreground is adopted. We caution 
that such estimates necessarily involve extrapolation from 
comparatively small numbers of BHBs. 

Could contaminating stars be responsible for the high 
velocity dispersion of 4.3±1.2 kms~^ reported bv lGeha et al.l 
(2009)? Figure [TT] shows the results of Monte Carlo simula- 
tions to gauge the importance of this contamination. The 
Sagittarius stream stars are assumed to have a velocity dis- 
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persion of 10 kms^^ and Segue 1 stars to have a velocity 
dispersion of 1 kms~^ (typical for a globular cluster). In 
the simulations, samples of 24 stars are generated, with the 
number of contaminating Sagittarius stream stars allowed 
to vary between and 10. The simulations are run for 10000 
iterations. The points on the plot are the average dispersion 
measured in the 10000 samples of 24 stars and the error bars 
represent the spread. 

The number of contaminating Sagittarius stars required 
to inflate the velocity dispersion from its true value (in the 
simulation) of 1 kms~^ to the value reported by Geha et 
al is surprisingly small - perhaps even one contaminating 
star is enough! Recall that we estimate the actual number 
of contaminants in the Geha et al sample is probably be- 
tween 4 and 8. Although the kinematic selection of Geha 
et al is sufficient to exclude stars from one of the wraps of 
the Sagittarius, it is not sufficient to exclude stars from the 
wrap moving with vqsr ~ 130 kms~^. A very modest con- 
tamination from this wrap can cause an anomalously high 
velocity dispersion to be reported. 



5.2 Tidal Influence 

Assuming that Segue 1 is on a Sagittarius-like orbit, its ra- 
dial and tangential velocities are sufficiently well constrained 
to ask what effect Galactic tides would have on the object. 
Figure [T2] shows the results of orbit integrations, in which 
a test particle is placed at the location of Segue 1 with a 
variety of initial velocities. The particle travels through the 
Milky Way potential as defined in Fellhauer et al. (2006) for 
roughly 2 periods. We then measure the perigalactic distance 
of the orbit as well as its eccentricity. The contours show 
the regions in which the points lie in the plane of perigalac- 
tic distance and eccentricity. We then calculate the tidal 
radiu s for a cluster of ma ss m travelling on such an orbit 
using llnnanen et al. I l|l983l ) 

= I ( (3T^) ^ 

Here, Rp is the perigalactic distance, e is the eccentric- 
ity, and Mp is the mass of the galaxy interior to Rp. In 
both panels of the figure, the radial velocity is Vr ~ 210 
kms~^ (heliocentric) or 118 kms~^ (Galactocentric). But, 
panel (a) shows the results with initial tangential veloci- 
ties corresponding to a typical Galactic range (75 kms~^ to 
225 kms~^), and with cluster masses of m = IOOOA/q and 
m= 10*^ Mq. Panel (b) assumes the tangential velocity is cho- 
sen to correspond approximately to that of the old leading 
arm 215 kms~^ Galactocentric). 

From both panels, we find that if Segue 1 is a globu- 
lar cluster with mass around IOOOA/q, then its tidal radius 
would be smaller than its observed half-light radius which 
implies that it could not survive on such an orbit for long 
and would be in the throes of destruction. If Segue 1 is a 
dwarf galaxy with mass of 10^ Mq, the estimated tidal ra- 
dius is much larger than the measured half-light radius and 
we would not expect to see any signs of tidal disruption 



6 DISCUSSION AND CONCLUSIONS 

The distinction between star clusters and dwarf galaxies is 
largely based on size. There is a factor of > 10 difference 
in the characteristic sizes of clusters and galaxies. Despite 
differences in the velocity dispersion, this translates into an 
order of magnitude difference in the estimated mass-to-light 
ratios, under the assumption of virial equilibrium. 

The example of Pal 5 shows the dangers of assuming 
virial equilibrium. Pal 5 has a half-light radius of ~ 20 pc, 
one of the largest of the Milky Way globular clusters. It is 
also well-known to be disintegrating under the Galactic tidal 
field, and so its half-light radius exceeds its tidal radius. The 
most distant stars are therefore already unbound and belong 
properly speaking to Pal 5's tidal tails. The assumption of 
virial equilibrium would lead to a mass-to-light ratio that is 
in serious error. As noted by Dehnen (2004), the final dis- 
ruption of Pal 5 happens very quickly, and most mass is lost 
during the final few percent of its lifetime. Another exam- 
ple of a dissolving Milky Way satellite is provided by Ursa 
Major II (UMa II). In the model of Fellhauer et al. (2007), 
the orbit of UMa 11 has pericentric and apocentric distance 
similar to that of the Sagittarius. The fate of UMa II is very 
similar to that of Pal 5, but the final disruption happens on 
a longer timescale, comparable to the orbital period of 1 Gyr 
(see Figure 6 of Fellhauer et al. 2007). Suppose we were to 
observe Pal 5 or UMa II at a future time equal to a fraction 
of their orbital period. Then, just a few faint stars will mark 
the nucleus, whilst any streams will become blurred and of 
lower contrast as they diffuse in the Milky Way potential. 
In other words, we will see something like Segue 1. 

Of course, to prove the hypothesis of extra-tidal stars 
for Segue 1, what is needed is wide-field kinem atical data. 
At th e moment, the only published study is bv lGeha et al.l 
(|2009l ). They measured the radial velocities of 24 stars and 
concluded that Segue 1 has a heliocentric velocity of ~ 206 
kms~^ and a velocity dispersion of 4.2 kms~^. The sys- 
temic velocity therefore coincides wi th both the predicted 
velocity of the Sagittarius stream (jFellhauer et al.l |2006| ) 
and the measured velocity as recorded in this paper. This 
not only strengthens the argument that Segue 1 is a star 
cluster, originating from the Sagittarius, but also raises the 
question of the levels of contamination in kinematically se- 
lected datasets around Segue I's location. In particular, such 
contaminants may be inflating the velocity dispersion, and 
hence the inferred mass-to-light ratio, of Segue 1. 

Taking advantage of the wide field of view of the 
Anglo- Australian Telescope/ AAOmega fiber-fed spectro- 
graph combination Q, Gilmore and collaborators have ob- 
served a 0.8° X 0.8° area around Segue 1. Observing pa- 
rameters were the same as those for the related study of 
Bootes 1 (see Norris et al. 2008). Targets are selected from 
the SDSS photometric data to be broadly consistent with 
the Segue 1 CMD published by Belokurov et al. (2007) and 
a radial velocity accuracy of < 10 kms~^ is attained. Analy- 
sis of this dataset is ongoing and will be presented elsewhere. 
Here, we are primarily concerned with existence of extratidal 
stars. The left panel of Figure [T51 shows the dataset plotted 
in the plane of heliocentric radial velocity and angular dis- 
tance from Segue 1. Applying a kinematic cut of 206 ± 20 

^ See |http: / /www.aao. gov.au/IocaI/www/ aaomega/ ] 
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Figure 12. Left: The results of orbit integration with initial velocity Vr ~ 210 kms"'^ and a spread of tangential velocities (75 kms~^ 
— 225 kms^^ and m = IOOOMq (left contour labels) and 10® Mq (right contour labels). The tidal radii implied by these orbits suggest 
that if Segue 1 was a globular cluster it should show signs of tidal disruption. Right: The same, but under the assumption that Segue 1 
is following the Sagittarius orbit of Fellhauer et al (2006). 
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Figure 13. Left: Heliocentric radial velocity of giant stars from the AAT survey plotted against angular distance from the center of 
Segue 1. The filled circles show 8 stars that have the same kinematics as Segue 1 (206 ± 20 kms-l). Right: The locations of the 8 stars 
are spread over a wide field of view. The grey contour lines show the density distribution around Segue 1 as taken from the middle panel 
of Figure |4] The black circle marks the half-light radius of Segue 1. The AAT candidates extend out to ~ 5 half-light radii. 



kms~ picks out the 8 likely members, which are reasonably 
well separated from the Galactic foreground populations. 
The right panel shows the location of the likely members 
as compared to the density distribution inferred from the 
optimal filter. Not only do the candidates extend well be- 
yond the half-light radius marked by the thin circle, but 
they also follow in broad outline the distribution suggested 
by the photometric analysis. This confirms the wide spa- 
tial distribution of stars with kinematics like Segue 1 in this 
field. 

Accordingly, we favour the scenario in which Segue 1 is 
a star cluster stripped early on from the Sagittarius galaxy. 
This provides a natural explanation for the fact that all four 
of the phase space coordinates that are measured coincide 
with those of the Sagittarius stream. What is needed to make 
the argument complete is a demonstration that the metallic- 



ity of S egue 1 is compatib le with an origin in the Sagittarius 
galaxy. ICeha et al.l (|2009r ) measured the metallicity of one 
star as [Fe/H] = -3.3, and the metallicity of a further 13 as 
[Fe/H] — -1.8. The latter is consistent with the metallicity 
of Sagittarius. The former is very surprising, as there are no 
Sagittarius clusters with a metallicity lower than -2.2. How- 
ever, the declination and radial velocity of Segue 1 place it 
in the very oldest wrap of material, and so it must have 
been stripped first and have resided in the very outskirts 
of the Sagittarius progenitor. This could explain variations 
in the metallicity of Segue 1 and other parts of the Sagit- 
tarius stream. Reg ardless of whether the metal-poor star of 
(|Geha et al.ll2009l ) is a member of Segue 1 or not, it signi- 
fies the detection of a substantial metallicity gradient in the 
Sagittarius progenitor, alth ough this is more extreme than 
has been measured so far (jChou et al.ll2007i ). In any case. 
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the claimed status of Segue 1 as the least luminous galaxy 
is very uncertain. 



ACKNOWLEDGMENTS 

We thank the referee, Rodrigo Ibata, for a careful reading of 
the paper. MNO is funded by the Gates Cambridge Trust, 
the Isaac Newton Studentship fund and the Science and 
Technology Facilities Council (STFC), whilst VB acknowl- 
edges financial support from the Royal Society. 

Funding for the SDSS and SDSS-II has been provided 
by the Alfred P. Sloan Foundation, the Participating In- 
stitutions, the National Science Foundation, the U.S. De- 
partment of Energy, the National Aeronautics and Space 
Administration, the Japanese Monbukagakusho, the Max 
Planck Society, and the Higher Education Funding Council 
for England. The SDSS Web Site is , http://www.sdss.org7] 



Tucker, D. L., et al. 2006, AN, 327, 821 

Willman, B., et al. 2005, ApJ, 629, L85 

Walsh S., Jerjen H., Willman B. 2007, ApJ, 629, L85 

Walsh S., Willman B., Sand D., Harris J., Seth A., Zaritsky 

D., Jerjen H. 2008, ApJ, 688, 245 
York D., et al. 2000, AJ, 120, 1579 



REFERENCES 

Adelman-McCarthy, J. K., et al. 2008, ApJS, 175, 297 

Belokurov, V., et al, 2006, ApJ, 642, L137 

Belokurov, V., et al. 2007, ApJ, 654, 2 

Carollo, D., et al. 2007, Nature, 450, 1020 

Chou, M.-Y., et al. 2007, ApJ, 670, 346 

Clem, J. L., 2005, PhD Thesis, University of Victoria 

Dehnen, W., Odenkirchen, M., Grebel, E. K., & Rix, H.-W. 

2004, AJ, 127, 2753 
Fellhauer, M., et al. 2006, ApJ, 651, 167 
Fellhauer, M., et al. 2007, MNRAS, 375, 1171 
Fukugita, M., Ichikawa, T., Gunn, J. E., Doi, M., Shi- 

masaku, K., & Schneider, D. P. 1996, AJ, 111, 1748 
Geha, M., Willman, B., Simon, J. D., Strigari, L. E., Kirby, 

E. N., Law, D. R., & Strader, J. 2009, ApJ, 692, 1464 
Gunn, J.E. et al. 1998, AJ, 116, 3040 
Gunn, J.E. et al. 2006, ApJ, 131, 2332 
Hogg, D.W., Finkbeiner, D.P., Schlegel, D.J., Gunn, J.E. 

2001, AJ, 122, 2129 
Ibata, R., Lewis, G. F., Irwin, M., Totten, E., & Quinn, T. 

2001, ApJ, 551, 294 
Innanen, K. A., Harris, W. E., & Webbink, R. F. 1983, AJ, 

88, 338 

Ivezic, Z. et al. 2004, AN, 325, 583 
Juric, M., et al. 2008, ApJ, 673, 864 

Martin, N. F., de Jong, J. T. A., & Rix, H.-W. 2008, ApJ, 
684, 1075 

Norris, J. E., Gilmore, G, Wyse, R. F. G, Wilkinson, M. I., 
Belokurov, V., Evans, N. W., & Zucker, D. B. 2008, ApJ, 
689, L113 

Odenkirchen, M., Grebel, E. K., Dehnen, W., Rix, H.-W., 

Cudworth, K. M. 2002, AJ, 124, 1497 
Odenkirchen, M., et al. 2003, AJ, 126 

Oke J. B., Gunn J. E. 1983, Ap J, 266, 713 

Padmanabhan, N., et al. 2007, arXiv:astro-ph/0703454" 
Pier, J.R., Munn, J. A., Hindsley, R.B, Hennessy, G.S., 

Kent, S.M., Lupton, R.H., Ivezic, Z. 2003, AJ, 125, 1559 
Schlegel, D.J., Finkbeiner, D.P., & Davis, M., 1998, ApJ, 

500 

Sirko, E., et al. 2004, AJ, 127, 899 
Smith, J. A., et al. 2002, AJ, 123, 2121 



© 2009 RAS, MNRAS 000.[HfT2l 



